Improving Transition Voltage Spectroscopy of Molecular Junctions 



Trocls Markusscn,^'^ Jingzhc Chen,^'0 and Kristian S. Thygcscn-'^ 

' Center for Atomic-scale Materials Design (CAMD), Department of Physics, 
Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark 
^Danish National Research Foundations Center of Individual Nanoparticle Functionality ( CINF), 
Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark 

(Dated: December 17, 2010) 

Transition voltage spectroscopy (TVS) is a promising spectroscopic tool for molecular junctions. 
The principles in TVS is to find the minimum on a Fowler-Nordheim plot where ln(//V^) is plotted 
against 1/V and relate the voltage at the minimum, Vmin, to the closest molecular level. Importantly, 
Knin, is approximately half the voltage required to see a peak in the dl/dV curve. Information about 
the molecular level position can thus be obtained at relatively low voltages. In this work we show 
that the molecular level position can be determined at even lower voltages, by finding the 

minimum of \n{I /V") with a < 2. On the basis of a simple Lorentzian transmission model we 
analyze theoretical ab initio as well as experimental I — V curves and show that the voltage required 
to determine the molecular levels can be reduced by ~ 30% as compared to conventional TVS. As 
for conventional TVS, the symmetry /asymmetry of the molecular junction needs to be taken into 
account in order to gain quantitative information. We show that the degree of asymmetry may be 
estimated from a plot of vs. a. 



PACS numbers: 73.63.Rt,73.23.-b 

I. INTRODUCTION 

The field of molecular electronics holds the promise 
of low-cost and flexible devices and potentially a con- 
tinued miniaturization of electronic circuits beyond the 
limits of standard silicon technologies The electronic 
properties of a molecular junction are to large extend 
governed by the position of the frontier orbitals with re- 
spect to the electrode Fermi levels, i.e. either the highest 
occupied molecular orbital (HOMO) or the lowest unoc- 
cupied molecular orbital (LUMO). One way of measuring 
the energetic position of the HOMO (or LUMO) level is 
from peaks in the AI /dV curve, where / is the current 
and V is the bias voltage. In the simplest picture, it re- 
quires a voltage olV ^ 2\Ef—£h\I&^o probe the HOMO 
position in a symmetric junction. However, in practice 
the junction will often suffer from breakdown before the 
peak in dl /AV is reached. 

Transition voltage spectroscopy (TVS) has recently 
been introduced as a spectroscopic tool in molecular elec- 
tronics^ and offers the possibility to probe the energy of 
molecular levels at relatively low voltages. Beebe et al^ 
showed that the HOMO position, determined from ultra- 
violet photo spectroscopy (UPS), scaled linearly with a 
characteristic voltage, Vmm obtained at the minimum of 
a Fowler-Nordheim graph, i.e. a plot of ln(//y^) against 
1/V. TVS is now becoming an increasingly popular tool 
in molecular electronics^ii"— . The great advantage of 
TVS is that the minimum voltage Knin ~ \Ef — £h\ in 
a symmetric junction is obtained at approximately half 
the voltage required for measuring a peak in dl/d^ i^^d^ 

The original interpretation of TVS is based on a tun- 
nel barrier modeU^. Within this model the transition 
voltage is reached when the tunnel barrier, due to the 
bias voltage, is shifted from a trapezoidal to a triangular 



shape. However, Huismann et alt^ pointed out that the 
barrier model in some respects is inconsistent with ex- 
periments and advocated that a Landauer transmission 
model with a single electronic level is more appropriate. 
In the Landauer model, the transmission function is as- 
sumed to be described by a single Lorentzian with a peak 
at the HOMO position. In a recent workii we studied 
a large number of molecular contacts by ah initio trans- 
port calculations, and showed that the simple Lorentzian 
transmission is a proper model, provided that the sym- 
metry/asymmetry of the junction is taken into account. 

When a bias voltage is applied across a molecular junc- 
tion, the two Fermi levels of the electrodes open a bias 
window in the gap between the HOMO and the LUMO. 
The dl/dV curve is peaked, when one of the Fermi lev- 
els is aligned with either the HOMO or the LUMO. At 
lower voltages, the Fermi levels go through the tail of the 
corresponding Lorentzian peak. The shape of the tail is 
determined by the level position, Eq, and the broadening, 
r, which is influenced by the contact geometries. As illus- 
trated in Ref. [ill the TVS minimum voltage is insensitive 
to the broadening provided that \eo — Ep\lV :s> 1. 

Our main focus in this paper is to discuss how to 
measure the HOMO or LUMO position at low voltages 
by exploiting characteristic features in the tail of the 
Lorentzian transmission peak. To this end we extend the 
TVS to a more general form, i.e., the objective function 
is switched to ln(//F") with a < 2, where we recover 
the usual TVS for a ~ 2. The minimum of In^/V") is 

denoted V^^. The properties of this function include: 

(i) the linear relationship between V|^",) and Eq still ex- 
ists which means any a value can in principle be used to 
predict Eq; (ii), a minimum value in the general TVS plot 
can always be obtained by tuning the a value; (iii), com- 
pared with conventional TVS, Vmin can be even lower 
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with smaller a. The last point is a substantial techni- 
cal improvement since the importance of TVS lies in the 
fact that the minimum in the Fowler-Nordheim plot is 
obtained at relatively low bias voltage. The analysis pre- 
sented in this article based on theoretical ab initio as 
well as experimental / — V^-curves shows that the volt- 
age required to determine the molecular levels can be 
reduced by 30%. As previously mentioned, the sym- 
metry/asymmetry of the molecular junction needs to be 
taken into account in order to gain quantitative infor- 
mation. We show that the degree of asymmetry may be 
estimated from a plot of V^l vs. a. 



II. MODEL 

As introduced in Ref. fill , the energy-dependent trans- 
mission coefBcient is assumed to be a Lorentzian 



T{E;eo,r) = 



f 



{E-eor + ryA' 



(1) 



where £o and F are the molecular level energy and broad- 
ening, respectively. / is a constant factor to account for 
multiple molecules in the junction, asymmetric coupling 
to the electrodes, etc. Our analysis will not be depen- 
dent on the actual value of /. We shall assume that the 
molecular level is the HOMO level, i.e. Eq < Ep. This 
is the typical case for thiol bonded molecules (as con- 
sidered in this work), where electron transfer from the 
metal to the sulfur end group shifts the molecular levels 
upward in energjsi^. At finite bias voltage, the molecu- 
lar level may be shifted relative to the zero-bias position. 
Such non-linear effects express themselves differently in 
symmetric and asymmetric junctions. We describe the 
degree of asymmetry by the parameter 77 e [—1/2; 1/2] 
such that the current is given by 
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T{E- £0 + VV,T) [fL{V/2) - fBi-V/2)] dE, 



(2) 

where fL/niV) = l/[exp{EF + eV)/kBT + 1] are the 
Fermi-Dirac distributions for the left and right contact, 
respectively. A symmetric junction has = 0, while a 
completely asymmetric junction has rj ~ ±1/2. 

Figure [1] (top) shows the I ~ V curve for a symmet- 
ric junction with a Lorentzian transmission function with 
£0 — Ep = —1 eV and F = 0.2 eV. The points marked by 
circles represent the minima of the \n{I/V") vs. 1/V 
plots shown in the bottom panel. The resonant trans- 
mission happens at the bias voltage V ~ 2Y , which 
can be reduced to 1.2 V by using the conventional TVS 
(a = 2.0). Using a smaller a value, the minimum can 
be obtained at even lower bias, i.e., 1.0 V and 0.6 V for 
a = 1.6 and a = 1.2 respectively. We thus see that by 
using a = 1.2, the voltage needed to measure the HOMO 
position is reduced by a factor of three compared to the 
resonant transmission and is only half the voltage re- 
quired using conventional TVS with a = 2. However, it 




Figure 1: Top: I — V curve from a symmetric junction and a 
Lorentzian transmission function with eo — Ep = — leV and 
r = 0.2 eV (sliown in tlie inset). Tlie three lower points on the 
curve mark the voltages, V^^^ at tlie the minima of Fowler- 
Nordheim plots \n{I /V") vs. 1/V. The numbers show the 
a values. The resonant transmission voltage is also marked. 
Bottom: Generalized Fowler-Nordheim plots where ln{I /V°') 
is plotted against 1/V. 



is also evident from Fig. [T] (bottom) that the minima in 
the TVS plots for small a values are more shallow and 
it might be difficult to accurately determine the mini- 
mum if the I — V data contains noise. Consequently, the 
uncertainty of TVS increases with decreasing a. 

We conclude the analysis of the Lorentzian transmis- 
sion model by considering the effects of asymmetry. Fig- 
ure [2] shows the ratio of the resonant level and minimum 
voltage vs. a obtained for a symmetric and completely 
asymmetric junction. For both junctions, the ratio in- 
creases by decreasing the a value. In accordance with 
Fig. [T]this implies that for a given value of £0, the gen- 
eralized TVS minimum, is reduced by decreasing a. 

From Fig. [5] it is seen that the reduction in V^l is the 
largest for a completely asymmetric junction (77 = 1/2), 
where the minimum voltage is reduced by a factor of 4 
in tuning a from 2 to 1.2. In the same a-range, the min- 
imum for a symmetric junction is reduced by a factor of 
2. 
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Figure 2: Ratio of the resonant level and minimum voltage, 



a obtained for a symmetric (77 = 0) and 



completely asymmetric (77 = 1/2) junction. 



III. AB INITIO CALCULATIONS 

Having established that the generalized TVS plot of 
ln(//y") can provide information of the molecular level 
position at significantly lower voltages than required by 
conventional TVS (using a = 2), we now turn to analyze 
more realistic I—V curves calculated from first principles. 

The current at finite bias voltage is calculated using 
DFT in combination with a non-equilibrium Green func- 
tion (NEGF) method as described in Ref. El. Our DFT- 
NEGF method is implemented in GPAW, which is a real 
space electronic structure code based on the projector 
augmented wave methodi^iii. We use the Perdew-Burke- 
Ernzerhof (FEE) exchange-correlation functionali^, and 
a 4 X 4 fc-point sampling in the surface plane. The elec- 
tronic wave functions are expanded in an atomic orbital 
basisii. In all calculations, the molecule and the closest 
Au layers are described by a double-zeta plus polariza- 
tion (dzp) basis set, while the remaining Au atoms are 
described by a single-zeta plus polarization (szp) basis. 

First, we investigate two proto- typical molecular junc- 
tions, Au-l,4-benzenedithiol-Au (Au-BDT-Au) and Au- 
phcnyl-STM (Au-Ph-STM) as shown in Fig. [3] (a) and 
(c) , representing the symmetric and asymmetric class re- 
spectively. From the transmission plot in Fig. [3] (b) we 
can see for the Au-BDT-Au junction, the HOMO level is 
about -1.2 eV and LUMO level is about 2.8 eV. Since the 
bias window opens symmetrically, the TVS measurement 
will be a measurement of the HOMO position. The trans- 
mission for the Au-Ph-STM junction, Fig. |3] (d), shows 
that the HOMO level is about -2.6 eV and LUMO level 
is about 2.8 eV. Since it is a highly asymmetric junction, 
the HOMO level will be more pinned to the Fermi level of 
left electrode, and when a positive bias voltage is applied 
the Fermi level of the right electrode will scan through 
the HOMO. The TVS plot will thus again be a measure- 
ment of the HOMO position. A detailed analysis of the 
transmission spectraii at finite bias voltages shows that 
the junction is not completely asymmetric, but rather 
has 77 = 0.33. This is due to the relative small vacuum 



gap between the right tip and the molecule causing a 
finite voltage drop over the molecule itself. Increasing 
the vacuum gap further would increase 77 but this would 
cause computational problems due to the finite range of 
the atomic orbital basis. 




Figure 3: Atomic structure of a symmetric Au-BDT-Au junc- 
tion (a) together with an asymmetric Au-Ph-STM junction 
(c). The corresponding transmission functions are shown in 
panel (b) and (d), respectively. 



From the calculated I — V curves (not shown) we find 
the generalized TVS minimum voltages, V^l^, for differ- 
ent a values. Based on curves like the ones shown in Fig. 
[2] for r) = 0.04 and 77 ~ 0.33 we obtain predicted values 
of the HOMO position at different a values. These pre- 
dicted HOMO values are shown in Fig. |3]vs. a and they 
are seen to converge to the "true" HOMO position as 
determined from the transmission plots (Fig. |3] (b) and 
(d)) as a is increased. For both junction, using a = 1.5 
is sufficient to predict the HOMO position with an error 
around 0.1 eV. This corresponds to a reduction in TVS 
minimum voltages by 23% and 31% for the symmetric 
(Au-BDT-Au) and asymmetric (Au-Ph-STM) junctions, 
respectively. 

It is well known that DFT is unable to accurately de- 
scribe energy gaps and level alignment of molecules at 
surface a^^i^" . The HOMO positions calculated in Fig. |4] 
can therefore not be expected to be in quantitative agree- 
ment with the experimental value. However, it is not our 
goal in this paper to use DFT to predict the HOMO 
position, but rather to show that the generalized TVS 
scheme can be applied on rcalisticly shaped I — V curves 
beyond simple models. In this respect, and since the ra- 
tio |£o — Ef\/V^I is independent on Eq, the results in 
Fig. |4] are still valid, irrispcctively of the inaccuracies in 
DFT. 

To get a broader overview of the feasibility of the 
generalized TVS we analyze a series of molecular junc- 
tions. All the original data we use here are the same as 
those in our previous paper—. Figure [5] shows the ratio 



4 




-Ph (TVS, 71=0.33) 
-BDT (TVS, 71=0.02) 



1.5 
a 



Figure 4: Predicted HOMO level of phenyl (Ph) and ben- 
zenedithiol (BDT) junctions in Fig. ^obtained from the TVS 
spectra using different a-values. The horizontal dashed lines 
mark the HOMO position as determined from the transmis- 
sion spectra in Fig. |3](b) and (d). 
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{sh — £^fI/v1 ?2 v s. asymmetry parameter rj. As also 
shown in Ref. Till , the ab initio data for a = 2 (panel 
a) agree well with the analytical model when the asym- 
metry is taken into account. Importantly, a similar good 
agreement between the DFT data and the model is seen 
for a = 1.5 (panel b). It thus seems to be of general valid- 
ity that quantitative information can be obtained using 
a = 1.5. Further decreasing a to 1.3 (panel c) leads 
to significantly larger deviations from the model calcula- 
tion. The deviation can be due to (i) the transmission 
peak differs from a Lorentzian; (ii) error in the asymme- 
try factor measurement because of non-linear effects; and 
(iii) for small values of a the determination of the TVS 
minimum voltage becomes uncertain. 

For practical purposes, it may be usefull to apply the 
approximate linear relations between \£H — Ep\/V^"J^ and 
77: 



\£h 
\£h 



EF\/Vit 

■ mill 



Ef\/V^ 



2.3 77 + 0.85 {a = 2) (3) 
3.7r/+l.l (a = 1.5). (4) 



If knowledge of the asymmetry parameter exists, the 
linear relations can be used to predict the HOMO (or 
LUMO) position using the general TVS scheme. We note 
in particular, that if the same predicted HOMO level is 
obtained from both Eq. ([3]) and (|4]) this is a strong indica- 
tion that the true transmission function is well described 
by a Lorentzian function and the above theory is appli- 
cable. On the contrary, if information about the HOMO 
position exist, e.g. from UPS experimental, Eqs. (I3|)-(|4]) 
can be used to determine the junction asymmetry. 

Based on the ab initio calculations, we conclude that 
a = 1.5 can be used in the general TVS measurement to 
obtain the HOMO (or LUMO) position without losing 
precision. Depending on the degree of junction asymme- 
try this reduces the required bias range with 20-30% as 
compared to conventional TVS. 
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Figure 5: Ratio of HOMO position and minimum voltage, 
\£h — Ef\/V^^I vs. asymmetry parameter rj. The minimum 
voltages are obtained with the generalized TVS scheme using 
a = 2.0 (a), a = 1.5 (b), and a = L3 (c). The ab initio 
calculations follow the Lorentzian model prediction for a — 
1.5 and a = 2.0, while substantial deviations are seen for 
a = 1.3. 



IV. ANALYSIS OF EXPERIMENTAL DATA 

In the previous section it was shown that the gener- 
alized TVS based on a Lorentzian transmission function 
provides a good description of ab initio calculations. In 
this section we proceed to analyze experimental I — V 
data using the generalized TVS. We obtain the experi- 
mental data from Beebe et al^ and Tan et al^. We man- 
ually read off the experimental data, and subsequently 
smoothcn the data by a Gaussion convolution. In Fig. 
[S] we show the obtained I — V data in panel (a) for a 
Au-anthracenethiol-Au junction^ and in panel (c) for a 
Au-terphenyl-thiol-Au junctionSi. Both experiments are 
performed with mono-thiol molecules using conducting 
atomic force microscope (AFM). The solid lines show the 
data after the smoothening. Panels (b) and (d) show the 
corresponding Fowler-Nordheim plot using a = 2. 

As shown in Ref. [ij and discussed above, quantita- 
tive use of TVS is only possible if knowledge about the 
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Figure 6: (a): I — V curve acquired from Ref . (circles) to- 
gether with a smooth curve obtained by convolution of the 
data points with a Gaussian function, (b) Fowler-Nordheim 
plot (using Q = 2) for the data points and the smooth func- 
tion, (c) and (d) show similar data obtained from Ref[^. 



1 

0.8 



0.6 



-^E 0.4 



0.2 




□ 



Tan data 
Beebe data 
-Model, T|=0 
Model, T|=0.25 
Model, Ti=0.5 



1.5 
a 



Figure 7: Ratio of V^'l/V^^^ obtained from the experimental 
data (symbols) and from the Lorentzian model (lines). The 
experimental data follow the model reasonably well for 77 = 
0.25. This asymmetry value is in good agreement with our 
previous calculationsii. 



junction asymmetry exist. However, it may be difficult 
to measure this quantity in experiments directly. While 
it is certainly evident that the I — V curves in Fig. [H] (a) 
and (c) are asymmetric, it is difficult to get a quantita- 
tive estimate of the asymmetry factor from these values 
alone. 

Within the Lorentzian transmission model V^,^",) is pro- 
portional to |£o — Ef \ provided that jso — Ep\ » F. The 
proportionality constant depends on a and on the asym- 
metry factor, 77. This implies that the ratio VJ^tn/^^n 
independent of both |eo — Ef\ and F. We may thus in 
principle determine ry from the dependence of V^'^\/V^^^ 

as function of a. In Fig. [7] we plot Kn?n/Knin 
for the two experimental data sets together with the 
model calculations based on a Lorentzian transmission, 
corresponding to three different asymmetry parameters. 
Wc first notice that the experimental data points fall 
within the model calculations corresponding to a com- 
pletely symmetric (77 = 0) junction and a completely 
asymmetric junction (r; ~ 0.5). This indicates that the 
true transmission functions indeed is well described by 
a Lorentzian shape. Secondly, we also see that both ex- 
perimental data sets follow approximately the model cal- 
culation corresponding to 77 = 0.25. We note that this 
asymmetry value is in good agreement with our previous 
calculationsii and we will use rj = 0.25 in the remaining 
analysis. 

From the generalized TVS minimum voltages ob- 
tained from the experimental data together with a model 
calculations like in Fig. [2] (for 77 = 0.25) we obtain pre- 
dicted HOMO values. The results are shown in Fig. [5] 
vs. a in panel (a) and vs. minimum voltage in panel 
(b). As previously seen for the ab initio calculations, the 
predicted HOMO values are reasonably well converged 
for a > 1.5 for both data sets. As seen in panel (b) this 
corresponds to minimum voltages reduced by ^-^30%. 
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Figure 8: Predicted values of Ieh — EfI vs. a (a) and vs. Vmin 
(b) for anthracene (Beebe data) and terphenyl (Tan data) 
junctions. The predicted EIOMO levels are reasonably con- 
stant for a > 1.5, thus indicating that the Lorentzian model 
is applicable in analyzing the experimental data. Panel (b) 
shows that using a ~ 1.5, the HOMO level can be determined 
at ~30% lower voltage than required for q = 2. 

From Fig. |8] we obtain approximate predicted val- 
ues of \eH — Ep\ « 0.85 eV and 1.0 eV for Beebe (an- 
thracene) and Tan (terphenyl) data, respectively. The 
HOMO positions of anthracene-thiol and terphenyl-thiol 
on a An surface were determined from ultraviolet pho- 
toelectron spectroscopy (UPS) to be Ep — er ^ 1.6 eV 
and Ep — £h = 1-7 cV, respcctivelyi^ The presence of the 
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AFM tip in the transport measurements will lead to a 
renormalization of the HOMO energy due to an image 
charge efFect^^ which is not present in the UPS measure- 
ments. Based on the method described in Ref. [21I we 
estimate the image charge interaction to be ^ 0.5 eV. 
Correcting the UPS values for eh by this value {eh = 
ejj + 0.5 eV) lead to Ep — Eh ~ 1.1 cV for the anthracene 
and Ep — eh ~ 1.2 eV for terphenyl, which agree with 
our TVS predicted HOMO values within 0.25 eV, corre- 
sponding to a deviation of ^ 25 %. The original interpre- 
tations of TVS assuming a one-to-one relation between 
the HOMO position and the TVS minimum voltage (i.e. 

\eh — Ef\ = V'n""^^) would have resulted in values of 
\eh — Ef\ = 0.62, 0.69 eV. These values differs signif- 
icantly from the image charge corrected UPS values by 
'-^ 75 %. The closer agreement between our predicted 
HOMO values and the UPS measurements signifies the 
importance of taking the asymmetry of the molecular 
junction into account. 

V. SUMMARY 

Transition voltage spectroscopy (TVS) is a valuable 
spectroscopic tool for molecular junctions as it offers the 
opportunity to probe molecular levels at relatively low 



voltages. In this work we have generalized the theory of 
TVS by showing that the minimum of ln(//y"), with 
a < 2, provides quantitative information of the HOMO 
position at lower bias voltages than required by conven- 
tional TVS, using a ~ 2. On the basis of a simple 
Lorcntzian transmission model we have analyzed theo- 
retical ah initio as well as experimental / — ^-curves and 
show that the voltage required to determine the molecu- 
lar levels can be reduced by ~ 30% as compared to con- 
ventional TVS. Importantly, the symmetry /asymmetry 
of the molecular junction needs to be taken into account 
in order to gain quantitative information. We believe 
the present work will be valuable for future theoretical 
as well as experimental work in the characterization of 
molecular junctions. 
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